Intrinsic Instability of Coronal Streamers 



Y. Chen 1 , X. Li 2 , H. Q. Song 1 , Q. Q. Shi 1 , S. W. Feng 1 , AND L. D. Xia 1 



ABSTRACT 

Plasma blobs are observed to be weak density enhancements as radially-stretched structures 
emerging from the cusps of quiescent coronal streamers. In this paper, it is suggested that the 
formation of blobs is a consequence of an intrinsic instability of coronal streamers occurring at a 
very localized region around the cusp. The evolutionary process of the instability, as revealed in 
our calculations, can be described as follows, (1) through the localized cusp region where the field 
is too weak to sustain the confinement, plasmas expand and stretch the closed field lines radially 
outwards as a result of the freezing-in effect of plasma-magnetic field coupling; the expansion 
brings a strong velocity gradient into the slow wind regime providing the free energy necessary 
for the onset of the subsequent magnetohydrodynamic instability; (2) the instability manifests 
itself mainly as mixed streaming sausage-kink modes, the former results in pinches of elongated 
magnetic loops to provoke reconnections at one or multi locations to form blobs. Then, the 
streamer system returns to the configuration with a lower cusp point, subject to another cycle of 
the streamer instability. Although the instability is intrinsic, it does not lead to the loss of the 
closed magnetic flux, neither does it affect the overall feature of a streamer. The main properties 
of the modelled blobs, including their size, velocity profiles, density contrasts, and even their 
daily occurrence rate are in line with available observations. 

Subject headings: Sun: corona — Sun: streamer — Sun: MHD instability 



1. INTRODUCTION 

Coronal streamers are quasi-stationary large- 
scale structures observed during solar eclipses 
or with coronagraphs, they are believed to be 
the outcome of a complex magnetohydrodynamic 
(MHD) coupling between the hot expanding coro- 
nal plasmas and the large-scale confining magnetic 
field (e.g., Pneuman & Kopp, 1971; Koutchmy & 
Livshits, 1992). Many streamers remain visible 
over a period of up to several months. Never- 
theless, various types of activities associated with 
coronal streamers are observed, including small- 
scale plasma blobs (Sheeley et al., 1997; Wang 
et al., 1998, 2000), streamer detachments with or 
without in/out pairs (Wang & Sheeley, 2006; Shee- 
ley & Wang; 2007; Sheeley, Warren & Wang 2007), 
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streamer blowout coronal mass ejections (CMEs) 
(Howard et al. 1985; Hundhausen 1993), and oth- 
ers (like streamer puffs (Bemporad et al. 2005), 
streamer deformations and deflections possibly in- 
duced by nearby coronal disturbances (e.g., Hund- 
hausen 1987; Sheeley, Hakala, & Wang 2000)). 

Among these activities, we are particularly in- 
terested in the formation of plasma blobs, which 
were discovered with the Large Angle Spectro- 
metric Coronagraph (LASCO) onboard the Solar 
and Heliospheric Observatory (SOHO) (Brueck- 
ner et al. 1995), first reported by Sheeley et al. 
(1997). They manifest themselves as weak density 
enhancements emerging from the elongated tips 
of coronal streamers and confined to the narrow 
plasma sheet without causing any permanent dis- 
ruption or dislocation of streamers. Since they are 
thought to flow passively along the bright plasma 
sheet and get dynamically coupled to the ambi- 
ent plasmas, blobs are regarded as tracers of the 
ambient flow velocity with only slight contribu- 
tion to the total slow wind mass flux. The esti- 
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mated density inside the blobs is weakly enhanced 
over its surroundings by less than, say, 10 %. The 
sizes, when blobs are first observed, are about 1 
R© in the radial direction and 0.1 R© in the trans- 
verse direction (corresponding to about 2° in lati- 
tudinal extension). The blob speeds generally in- 
crease from - 250 km s _1 when first observed 
with LASCO C2 in the region 2 - 6 R Q (solar 
radii) to about 250 - 450 km s _1 in the region 
20 -30 R Q (Sheeley et al. 1997). It was shown 
by Wang et al. (1998) that the occurrence rate 
of blobs is about four per day along an edge-on 
streamer in the April of 1997, the quiet period of 
solar activities. During times of solar eruptions, 
the number, size and brightness of blobs generally 
increase (Wang et al. 2000). 

Theoretical modelling of coronal streamers 
started early in the 1970s (e.g., Pneuman & Kopp, 
1971). The classical model developed by Pneuman 
& Kopp treated the streamer as a magnetostatic 
structure in a bipolar configuration with a thin ax- 
isymmetric current sheet. The sheet extends from 
the tip of the closed streamer arcades to infinity, 
which separates the open field lines with opposite 
directions. Subsequent models presented similar 
morphological picture of streamers, many of them 
concentrate on the heating and acceleration of 
the corona and solar wind plasmas emphasizing 
the steady or quasi-steady feature of the streamer 
structure (e.g., Steinolfson et al., 1982; Wang et. 
al., 1993, 1998; Suess et al., 1999; Chen & Hu, 
2001; Hu et al., 2003a; Suess & Nerney, 2006). 
Some authors did mention that it may be difficult 
to obtain an exact fully- converged steady state 
of the coronal streamer. For instance, Washimi 
et al. (1987) found that the modelled values of 
the velocity and density have a consistent tem- 
poral variation of about 1 % in magnitude per 
hour; Suess et al. (1996) demonstrated that the 
streamer is subject to slow and continuing evapo- 
ration till all the closed field lines are fully opened 
if heated with a volumetric heat source. Endeve 
et al. (2003, 2004), on the other hand, showed 
that coronal streamers can undergo eruptions if a 
specific heat source deposited inside the closed ar- 
cades is added to protons rather than electrons. It 
should be pointed out that the numerical resolu- 
tion along the latitudinal direction in their studies 
is taken to be 2.77° between adjacent grid points, 
which is even larger than the latitudinal size of 



the blobs observed by Sheeley et al, as introduced 
previously. Therefore, in these models it is impos- 
sible to resolve the blob structures as well as their 
formation process. 

In Einaudi et al. (1999), the formation of blobs 
is simulated as magnetic islands resulting from the 
nonlinear development of the tearing mode insta- 
bility along a hypothetical one dimensional Har- 
ris neutral sheet. The sheet is assumed to be 
embedded in a prescribed fluid wake, as a sim- 
plification of the realistic configuration associated 
with the streamer belt above the cusp. Refine- 
ments of this oversimplified incompressible model 
considered the effects of compressibility and in- 
homogeneous density distribution in relevance to 
the dense plasma sheet (Einaudi et al., 2001). By 
further considering the two-dimensional configu- 
ration with open and closed magnetic field lines 
and the specific cusp geometry, Lapenta & Knoll 
(2005) concluded that the reconnections forming 
the plasma blobs are driven by converging flows 
at the cusp region with reconnection rate mainly 
determined by the driving flows instead of resis- 
tivity, allowing their results to be extended to 
realistic coronal Lundquist numbers. However, 
their model is based on the Cartesian geometry 
with prescribed wake velocity field, and the im- 
posed streamer morphology is not result from the 
self consistent dynamical equilibrium between the 
plasma expansion and magnetic confinement. In 
this paper, we present our theoretical endeavor 
in improving these relevant previous models for a 
better understanding of the formation of streamer 
blobs. The improvements are three-folds. Firstly, 
our calculations are conducted in the spherical ge- 
ometry, which is somewhat more realistic com- 
pared with the planar cartesian geometry. Sec- 
ondly, the streamer system in our model is self- 
consistently determined by the MHD coupling be- 
tween plasmas and magnetic fields. Lastly, we 
take advantage of a much higher numerical res- 
olution to better resolve the physical processes in 
the blob forming region. We show that the blobs 
are the consequence of an intrinsic instability of 
coronal streamers associated directly with the lo- 
calized cusp region. In the following section we 
reveal the details of our theoretical model with 
initial and boundary conditions. The obtained so- 
lution for the streamer solar wind-blob system is 
shown in section 3. A brief summary of this article 
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with discussion is provided in the last part of this 
article. 

2. Theoretical models 

Assuming the corona to be axisymmetric as 
in most previous studies, we solve the full- 
set two-dimensional and three-component (2.5(1) 
ideal MHD equations in the spherical coordi- 
nate system (r, #,</?). A magnetic flux function 
ip(t, r, 0) is introduced to express the magnetic 
field B(B r ,B e ,B (p ) as 

The deduced equations are given in Hu (2004) and 
are not to be repeated here. The equations are 
solved numerically in a domain of i? <r< 30Rq 
and < 9 < 90°, which is discretized into 300 x 
632 grid points. Note that we also assume the 
system to be symmetric about the equator, thus 
the parameters in the other half of the latitudi- 
nal domain is derived from the obtained solution. 
The grid spacing increases according to a geomet- 
ric series of a common ratio 1.015 along the radial 
direction from 0.005 Rq at the solar surface to 
0.437 Rq at the top boundary. In the latitudinal 
direction, we first adopt a uniform mesh between 
9 = 90° to 9 — 30° with grid spacing taken to 
be 0.1°, then we let the grid spacing increase ac- 
cording to a geometric series of a common ratio of 
1.134 from 9 — 30° to the pole. Such arrangement 
allows us to resolve the thin current sheets above 
the streamer cusps and the detailed dynamics of 
plasma blobs with high resolution. The multistep 
implicit scheme (MIS) developed by Hu (1989) is 
used to solve the MHD equations. The scheme has 
a relatively low numerical dissipation, which has 
been employed to deal with various issues in space 
physics, like the two-dimensional corona and solar 
wind solution driven by Alfvehic turbulence (Chen 
& Hu, 2001; Hu et al., 2003a, Li et al. 2004, 2006 ), 
and the flux rope catastrophe model for CMEs (Hu 
et al., 2003b; Chen et al., 2006, 2007). Since we 
are solving the equations under the framework of 
ideal MHD, any reconnections involved in the nu- 
merical solutions stem from numerical rather than 
physical resistivity. Therefore, the low dissipative 
feature of the adopted numerical scheme is critical 
to the issues we are concerned. 



The heating and acceleration mechanisms of the 
corona and solar wind plasmas, which still remain 
unresolved at the present time, are approximated 
by using a polytropic process with the polytropic 
index taken to be 7 = 1.05. To start our mod- 
elling, the initial magnetic field is taken to be a 
quadrupolar one with the azimuthal component 
set to be zero and the magnetic flux function given 
by Antiochos et al. (1999) as follows, 

^ , r sin 2 0R Q (3 + 5 cos 20) sin 2 OR^ 3 
^(r, 9) = M ^ + ^ 



2r 3 



(2) 

where ip c = , the unit of magnetic field 

strength Bo = 0.833 Gauss, the unit of magnetic 
flux function ipo = 4.04 x 10 13 Wb, and the mag- 
netic field strength at the pole B p is taken to be 8 
Gauss. It is known that a neutral point exists in 
this potential field at rjy = V^Rq along the equa- 
tor with ip = 0.3697^0- The reasons we adopt 
this specific topology, rather than the often used 
bipolar topology for this study are given as fol- 
lows. First, the triple-streamer system is more 
realistic, especially when one models the corona 
at high solar activity. As a matter of fact, even 
at solar minimum, the corona is often observed 
to have a few rather than one-single well-defined 
streamers. Second, this allows us, in future de- 
velopment of the model, to extend the break-out 
CME model originally proposed by Antiochos et 
al. (1999) to the case including the effect of the so- 
lar wind. Further elaboration concerning the rel- 
evance to the CME modelling will be given in the 
last section. Third, by focusing on the dynamics of 
the side streamer which is fully located inside the 
calculation domain, we avoid the influence of the 
symmetric boundary conditions along the equator. 

In Figure la we show the magnetic field lines 
of the initial quadrupolar field from 1 to 10 it!©, 
represented by the contours of the magnetic flux 
function with the contour level separated by ip = 
0.05^0- On the solar surface, the magnetic flux 
function first increases from at the pole to 1.2 
ipo near the center of the side arcade, and then 
decreases again to at the equator. The sepa- 
ratrix in the initial quadrupolar field specified by 
ip = 0.3697^0 is shown as red curves. It is easy to 
see that the initial magnetic field consists of four 
flux system including the central arcade system 
straddling the equator, two side ones centered at 
9 = 45° and 135°, and an overlying polar flux sys- 
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tern. As will be clearly illustrated in the following 
section, the initial magnetic topology is subject 
to significant change by the plasma-field coupling 
process. 

At the lower boundary (r = li?©), the mag- 
netic flux function ip is fixed to the values given 
by the above equation, the plasma temperature is 
set to be constant at 2 MK. The density, which is 
assumed to be dependent on the latitude, is 10 8 
cm -3 and 1.5xl0 8 cm -3 in the region = [0, 30°] 
and [60°, 90°], respectively, and increases linearly 
from = 30° to 60°. Note that although the as- 
sumption of an inhomogeneous density distribu- 
tion at the lower boundary results in a discontinu- 
ous latitudinal pressure gradient there, essentially 
similar physical processes occur when assuming 
a homogeneous density distribution, as revealed 
by calculations not presented here. The radial 
and latitudinal velocity components at the lower 
boundary are obtained by the mass conservation 
and the parallel condition given by V \\ B. The az- 
imuthal components of the velocity and magnetic 
field are set to be zero. All parameters are extrap- 
olated to get the conditions at the top boundary, 
and symmetric conditions are employed at both 
the equator and the pole. We are aware that 
the boundary conditions are somewhat arbitrarily 
chosen in order to obtain a solution with reason- 
able streamer structures and solar wind param- 
eters, they are still in line with available obser- 
vational constraints. The initial solar wind con- 
dition is an educated guess which does not have 
impacts on the final large scale converged numeri- 
cal solution as confirmed by the calculations. For 
example, the radial velocity is set to be a function 
of radial distance increasing monotonically from 
1.1 km/s at the bottom to 400 km/s at the top. 
The contour map of the two-dimensional distribu- 
tion of the prescribed radial velocity field is de- 
picted in Figure la. With the initial and bound- 
ary conditions described above, the coupling be- 
tween the plasma dynamical motion and the coro- 
nal magnetic field is self-consistently taken care 
of using our numerical code. The calculation fi- 
nally converges to a quasi-steady solution featured 
by a triple-streamer, slow solar winds emanating 
from the poles and between the closed arcades, 
and plasma blobs continuously released from the 
streamer tips, as will be depicted in the following 
section. 



3. The triple-streamer solar wind-blob 
system 

To save computation time required for the so- 
lution to converge, the calculations were first con- 
ducted in the same domain but with much less grid 
points in the latitudinal direction than described 
in the last section till a quasi-steady state was 
reached, which was redistributed to the present 
denser grid points. In the case with lower resolu- 
tion the number of grid points in the latitudinal 
direction is taken to be 150 with grid spacing uni- 
formly distributed between = 90° to = 30° 
with a resolution of 0.5°, and increasing according 
to a geometric series of a common ratio of 1.046 
from = 30° to the pole. The number and dis- 
tribution of grid points in the radial direction are 
taken to be the same in all calculations. The ob- 
tained solution was then allowed to evolve self- 
consistently for a long time, e.g., for a few hun- 
dred hours. The magnetic field lines with contour 
map of the radial velocity field of a solution at 
an arbitrarily-selected instant are shown in Figure 
lb. For the purpose of illustration, the time of this 
solution is set to be zero (i.e., t=0). As mentioned 
previously and clearly indicated by Figures la and 
lb, the original quadrupolar topology of the mag- 
netic field is reshaped significantly by the presence 
of the solar wind plasmas. The most apparent 
change is the formation of the triple-streamer sys- 
tem with solar wind plasmas streaming into inter- 
planetary space from both polar regions and the 
narrow regions between the central and the side 
streamers. In the outer part of the original side 
arcade system a few field lines are opened by the 
solar wind plasmas, while a few other field lines 
reconnect at the current sheet along the equator 
and close back resulting in an increase of the closed 
magnetic flux of the central arcades. This causes 
the disappearance of the original neutral line at 
the intersection of the red curves (ip = 0. 3697^0 ) 
in Figure la, and has a profound impact on the 
break-out model for CMEs proposed by Antiochos 
et al. (1999) based on the field topology shown in 
Figure la. We will return to this point in the last 
section. The field lines with the same value of ip 
are also plotted in Figure lb to Id for comparisons. 

In Figures lc and Id we illustrate the mag- 
netic field topology and the radial velocity field 
at t=21.5 and 30 hours, respectively, so as to in- 
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dicate the long-term evolution of the dynamical 
streamers. To show features in the contour map 
of the radial velocity, the magnetic field lines are 
omitted in one of the hemispheres in Figure lb 
to Id. An accompanying mpeg animation of the 
temporal evolution of the solution from t=0 to 48 
hours with frames separated by 10 minutes in time 
interval is provided with the online version of this 
paper. The solution at t=21.5 hours will be re- 
shown in the following figure to reveal the streamer 
dynamics in more detail. Note that the equator is 
taken as a boundary with symmetric conditions, 
we are aware that such treatment certainly has 
impacts on the dynamics associated with the cen- 
tral streamer. We will revisit this issue again as 
we proceed. Therefore, in this article we concen- 
trate only on the dynamics atop the side streamer. 
Yet, the modelling results in the central streamer 
region can still be used for the purpose of compar- 
ison. 

It is obvious from Figure 1 and the online an- 
imation that the streamer is in a constant state 
of motion with small magnetic/plasma islands re- 
leased continuously from the elongated tips of the 
side streamers, despite the long-term large-scale 
stability of the overall morphology. As indicated 
by the solutions, generally speaking, the islands 
have no regular shapes, the latitudinal width is 
about 0.1 to 0.2 i? , and the radial length about 
1 to 1.5 Rq when the islands are first separated 
from the main streamer. The plasma inside the 
islands is denser than the surrounding solar wind 
plasmas. For example, for the solution at t=23.5 
hours, which will be further illustrated in the fol- 
lowing figures, the center of an island just released 
is located at about 3.6 Rq. The density at this 
distance is about 1.5 xlO 6 cm -3 inside the island, 
and 1.4 xlO 6 cm -3 in the nearby equatorward so- 
lar wind and about 1.0 xlO 6 cm -3 in the nearby 
polarward wind. The quantities of our modelled 
magnetic islands are basically consistent with the 
LAS CO observations for streamer blobs reported 
by Sheeley at al. (1997) and Wang et al. (1998), 
and summarized in the introduction part of this 
article. We therefore regard them as streamer 
blobs hereinafter. The white field line delineating 
the outer boundary of streamer blobs in Figure 1 
and the online animation is given by ip = 0.65^o- 
To reveal more details of the blob dynamics, we 
also plot three more adjacent magnetic field lines 



with = 0.655,0.66,0.665^0 m green, red, and 
yellow colors, referred to as line A, B, and C, re- 
spectively. It can be seen that among them line 
C shows no apparent motion during the process, 
indicating that the blob formation involves only 
the localized cusp and current sheet region. The 
process does not result in the loss of the closed 
magnetic flux of the side streamer. 

As seen from the online animation, daughter 
blobs are found to form as a result of the breakup 
of a longer mother blob. After the formation, 
blobs change their length and shape continuously 
during the outward propagation. They usually get 
wider as a result of the relaxation of magnetic 
tension of reconnected field lines. It is also ob- 
served that two approaching blobs with differential 
speeds may get merged into one blob structure at 
certain distances. After checking the animated so- 
lution carefully, we get an average occurrence rate 
of about 6-8 blobs per day. Note that this estimate 
has excluded the blobs that are possibly too small 
to be measurable with current coronagraphs. The 
propagation velocity profiles of blobs, as estimated 
from our numerical results, are found to basically 
follow that of the surrounding solar wind plasmas. 
These modelling features can be regarded as in a 
remarkable agreement with available observations 
at the quiet time of solar activity considering the 
simplicity of our model. 

As clearly seen from the velocity maps depicted 
in Figure lb to Id, there exists apparent velocity 
shear on the two sides of the blob forming region. 
It is natural to surmise that the formation of blobs 
is related to the presence of the strong velocity 
shears that are capable of driving various modes 
of streaming MHD instabilities. To get deeper in- 
sights of how blobs are formed, in Figure 2 we 
show 10 image frames spanning from t=19 to t= 
23.5 hours, separated by 30 minutes in time inter- 
val with time increasing from left to right and top 
to bottom. All images are the side streamer part 
cropped from their corresponding global versions 
which are similar as that shown in Figure 1. For 
example, Figure 2f is corresponding to Figure lc 
(t=21.5 hours). From Figure 2a, we see that a nar- 
row long blob (referred to as blob I) is just formed 
from magnetic reconnection and starts its outward 
propagation. About 4 to 4.5 hours later, another 
blob is released from the same streamer (referred 
to as blob II) following the trend of its predecessor. 
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Therefore, the images collected in Figure 2 delin- 
eate a complete cycle of blob formation of typical 
configuration. Below we shall describe the evolu- 
tionary process and physical interpretation of the 
formation of blob II by illustrations with Figures 
2 and 3. 

For a visual guide, in Figure 2 we add three 
black arcs centered at the sun with radius being 
2.35, 2.6, and 3.6 Rq. We see from Figure 2a that 
the streamer cusp is reformed at a lower distance 
of about 2.6 Rq upon the disconnection of blob 
I. The magnetic field at the newly reformed cusp 
region is apparently too weak to sustain the con- 
finement of the plasmas at coronal temperature, 
the plasmas expand and stretch the freezing-in 
closed field lines radially outwards. Consequently 
the streamer tip elongates as clearly seen from the 
upper panels of Figure 2. After a continuing ex- 
pansion for about 2 to 2.5 hours line A (in white 
color with ip — 0.655^o) starts to pinch at its mid- 
dle part. It is apparent that the physics accounting 
for the pinching of field lines drives the following 
reconnection shown in Figure 2i and 2g. To check 
the physical conditions at the initial stage of the 
line pinch, in the left panel of Figure 3 we plot the 
radial components of velocity and magnetic field 
vectors (v r and B r ) together with the Alfven speed 
(va) for the solution shown in Figure 2f (t=21.5 
hours) as a function of co-latitude (A = 90° — 0). 
The profiles are plotted at the 3 heliocentric dis- 
tances used for the visual guide, i.e., 2.35 (solid), 
2.6 (dotted), and 3.6 (dashed) Rq. 

From the v r profile, we see that there exists a 
strong velocity jump between the surrounding so- 
lar wind plasmas along open field lines and that 
inside the elongated closed loops reaching up to 
nearly 150 km s" 1 at 2 to 3 solar radii. In the 
shear region the magnetic field changes its direc- 
tion, and the overall strength is relatively weak 
compared to that in the surrounding solar wind 
regimes. Note that B r at 2.6 Rq (the dotted line) 
changes its polarity 3 times in the shear region. 
This feature is an aftermath of blob I, and will be 
interpreted later together with the physical ori- 
gin of the small concave morphology of lines A 
and B at the streamer tips upon the formation of 
blobs I and II. The Alfven speed reaches the mini- 
mum value inside the velocity shear, as a result of 
the weak magnetic field there. Generally speak- 
ing, from 2.35 to 3.6 Rq the Alfven speed is less 



than 75 km s _1 in the solar wind flowing between 
the central and the side streamers, much higher in 
the polar wind regime, and less than 50 km s _1 
inside the velocity shear region. 

Velocity shears in plasmas can generate Kelvin- 
Helmholtz (KH) instabilities (e.g. Miura & Pritch- 
ett 1982). The physical configuration as shown in 
Figure 3a to 3c is somewhat similar to that stud- 
ied previously by, e.g., Lee, Wang & Wei (1988) for 
streaming instabilities (more widely known as KH 
instability), although they investigated the situ- 
ation with velocity maximum located inside the 
current sheet region, on the contrary to our situ- 
ation. Such discrepancy can be easily reconciled 
if one applies a coordinate transformation from 
their configuration into a moving reference sys- 
tem with an appropriate velocity, e.g., the max- 
imum velocity in their configuration. With the 
transformation, the real part of the frequency is 
doppler-shifted, and the instability criteria, eigen- 
mode profiles, and growth rate remain unaffected. 
Using their linear model of streaming instability, 
Lee et al. (1988) determined the eigenmode pro- 
files and critical conditions (see Equations (41) 
and (42) in their article) for both the sausage and 
the kink modes. To take advantage of the deduced 
conditions, we need to specify the Alfven speed in 
the shear region and in the surrounding plasmas, 
the number density, the wave number, and the half 
width of the shear region. Referring to the param- 
eter profiles plotted in Figure 3a to 3c, we assume 
an Alfven speed to be 20 km s _1 inside and 75 
km s _1 outside the shear region, the half width of 
the shear to be 0.05 Rq, and the density to be dis- 
tributed homogeneously. After simple calculations 
with Lee et al.'s critical conditions, we find that 
in the hypothetical configuration for the specified 
wave length of 1 Rq (typical length of a blob), the 
velocity gradient should be no less than 158 km 
s _1 for the sausage mode and 100 km s _1 for the 
kink mode to take place. By checking Figure 3a, 
we find that the total velocity jump presented in 
our solution at the time field line A starts to pinch 
can indeed reach up to 150 to 160 km s _1 . This 
implies that both the sausage and the kink modes 
of streaming instabilities can occur in our config- 
uration. We emphasize that the above estimate 
is rather rough and carried out with a much sim- 
pler configuration than that represented by Fig- 
ure 3a to 3c. It is well known that the sausage 
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mode results in the variation of the current sheet 
thickness to form a sausagelike configuration, and 
the kink mode produces a wavy structure keeping 
the current sheet thickness relatively constant. By 
checking the online animation again, it is appar- 
ent that both the sausagelike and wavy structures 
are present in the blob forming region. Thus, the 
occurrence of the mixed sausage-kink mode insta- 
bility in our solutions is actually self-evident. On 
the other hand, the symmetric boundary condi- 
tions adopted at the equator effectively prohibit 
the growth of the kink mode. Hence only stream- 
ing sausage mode can be observed atop the central 
streamer. 

From the above text, we propose that the pinch 
of field lines is caused by the streaming sausage 
mode instability. The instability, ideal in nature, 
is unable to produce reconnect ions. Yet, with the 
nonlinear growth of the instability, the distance 
decreases and the current density increases con- 
sistently in between the elongating field lines till 
whatever dissipation mechanism sets in to cause 
reconnection. Therefore, we suggest that the re- 
connections accounting for the blob separations 
are driven by the nonlinear development of the 
sausage mode. The reconnections could be caused 
by the well-known resistive tearing mode instabil- 
ity, or simple magnetic diffusion process across the 
narrow layer with enhanced current density. To 
determine the exact process acting in this highly- 
nonlinear MHD simulation, one needs to consider 
a more sophisticated and physical model of resis- 
tivity and a finer numerical resolution to reduce 
the numerical resistivity as low as possible. This 
is generally rather robust and computationally ex- 
pensive, and remains unresolved in the current 
study. As apparently seen from the multi-blob 
formation in our numerical solution, reconnections 
can take place simultaneously at various locations. 
The wavelength of the sausage mode is reflected 
by the length of the blob just formed, which is 
about 1.0 Rq to 1.5 Rq. The upper limit of the 
growth time of the sausage mode can be estimated 
by measuring the time span from a field line starts 
to pinch to a blob is finally released, which is about 
2 to 3 hours according to Figure 2e to 2j. 

The sausage mode instability is energized by 
the free energy associated with the velocity shear 
present in the blob forming region. To see how ve- 
locity shear evolves after the formation of a blob, 



in the right panels of Figure 3 (3d - 3f) we plot the 
latitudinal profiles for the solution shown in Fig- 
ure 2j (t=23.5 hours) at the same three distances 
as used in the left panels. From Figure 2j, it can 
be seen that the break point of the streamer is lo- 
cated at about 2.6 i? , the profiles at this distance 
are plotted as dotted lines. The velocity shear al- 
most vanishes crossing the streamer break point, 
as the aftermath of the sausage-mode instability. 
It can be seen from the solid lines plotted in both 
Figures 3a and 3d, strong velocity shear is always 
present at a lower distance of 2.35 Rq. Such ve- 
locity shear is along the streamer leg between the 
open and closed magnetic field lines, where the lo- 
cal magnetic field is relatively strong and a KH 
instability is unable to develop. In the following 2 
to 3 hours of continuous streamer elongation after 
t=23.5 hours, the strong velocity shear is brought 
into the slow wind regime again where the field 
strength gets lower so that the critical conditions 
for streaming instabilities can be satisfied. 

Similar to the profile for 2.6 Rq plotted in Fig- 
ure 3b, the radial magnetic field at about 2.35 Rq 
in Figure 3e also presents a triple change of po- 
larity. This feature can be straightforwardly ex- 
plained by the concave- inward morphology of line 
B and its adjacent lines not plotted here. The for- 
mation of the concave-inward morphology can be 
interpreted as follows. The causes are different for 
different field lines. For the field line that is just 
broken by reconnections, like line B (in red color) 
in Figure 2i, the morphology is a result of the 
relaxation of intensified magnetic tension of the 
over-stretched field lines upon reconnections. For 
the field line overlying the lines just reconnected, 
like line A in Figures 2a and 2i, the morphology is 
a result of magnetic pressure gradient force which 
points inwards as a result of the rapid shrinking 
of those lines just reconnected (like line B). 

As a summary of this section, the formation 
of streamer blobs can be regarded as consisting 
of two successive processes, one is the plasma- 
magnetic field expansion through the localized 
cusp region where the field is too weak to sus- 
tain plasma confinement, the other is the onset 
of the sausage mode instability that finally trig- 
gers magnetic reconnections producing separated 
plasma blobs. The former process brings strong 
velocity shear into the slow wind regime, and pro- 
vides required free energy for the following stream- 
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ing instabilities. By checking the time frames col- 
lected in Figure 2 for a complete cycle of the for- 
mation of a typical blob, we find that generally the 
streaming sausage mode instability sets in after 
the expansion continues for about 2 to 2.5 hours. 
The expansion continues till the nonlinear devel- 
opment of the instability finally breaks the elon- 
gated loops, which takes another 2 to 3 hours, 
producing a streamer blob and a 'new' coronal 
streamer with a lower cusp. The whole process 
originates from the peculiar magnetic topological 
feature at the cusp region intrinsically associated 
with a coronal streamer. We therefore term the 
described process as intrinsic instability of coronal 
streamers. Although the instability is intrinsic, it 
does not lead to the loss of the closed magnetic 
flux, so it does not influence the overall feature of 
a streamer. 

4. Conclusions and discussion 

It is well known that a coronal streamer sys- 
tem is resulted from the dynamical equilibrium 
between plasma expansion and magnetic confine- 
ment. Despite the overall long-term stability of 
the morphology, it is suggested that the streamer 
is subject to an intrinsic instability originating 
from the peculiar magnetic topological feature at 
the cusp region. The whole process of the insta- 
bility consists of two successive processes, one is 
the plasma-magnetic field expansion through the 
localized cusp region where the field is too weak to 
maintain plasma confinement, the continuing ex- 
pansion brings strong velocity shear into the slow 
wind regime providing free energy necessary for 
the onset of a streaming sausage mode instability 
which causes pinches of magnetic field lines and 
drives reconnections at the pinching points to form 
separated magnetic blobs. 

Although our recent two-dimensional models 
for the corona and solar wind have considered the 
role of Alfven waves in the formation of the bipolar 
coronal streamer and the solar wind as a two-fluid 
and two-state plasmas (Chen & Hu, 2001; Hu et 
al., 2003a), in this paper we still take advantage 
of the polytropic assumption of the heat transport 
process for simplicity. These previous streamer 
and solar wind models were aimed at investigat- 
ing the heating and acceleration of the corona 
and solar wind plasmas emphasizing the large- 



scale stability of the corona streamers rather than 
the detailed dynamics of blobs. To do this, a cer- 
tain amount of numerical dissipation was applied 
to stabilize the calculations. As a drawback of 
this treatment, the dynamical processes associated 
with blob formations have been largely smoothed 
out in these previous models. Note that although 
no explicit heating is applied to the plasma in the 
present model, a polytropic index of 7 = 1.05 and 
a temperature of 2 MK at the inner boundary re- 
sults in the high-temperature coronal plasma with 
a nature to expand. In a future extension to the 
present model, one may let 7 = | and employ a 
more realistic energy balance equation to model 
the plasma heating process. This will change the 
details of the plasma-magnetic coupling process, 
and thus the streamer morphology as well as the 
solar wind and blob properties. However, we be- 
lieve that the overall physics of blob formation sug- 
gested above, namely, the instabilities are essen- 
tially due to the presence of an intrinsic streamer 
cusp with the dynamics determined by the elon- 
gation of the freezing-in plasma-magnetic arcades, 
which brings the relatively slow confined plasmas 
into the solar wind regime resulting in enough ve- 
locity shear to trigger the streaming sausage mode 
instability, remains largely unaffected. 

In the study by Lapenta & Knoll (2005), blobs 
are interpreted as a result of reconnections along 
open current sheet structure beyond the streamer 
cusps. This can be seen by comparing their Fig- 
ure 6a and Figure 6d (also 7a and 7c, 7d and 
7f). At the end of their case calculation there is 
an apparent increase of the closed-flux inside the 
streamer. The increase can be more than 50-100% 
of the original total closed flux contained by their 
streamer if the magnetic flux function along ad- 
jacent field lines plotted in their figures is equally 
separated. Thus, the streamer itself will grow sig- 
nificantly and continuously after blobs form one by 
one. Although the magnetic field is not directly 
observable in the corona, this shall be inconsis- 
tent with observations which indicate that blobs 
have only minor contributions to the solar wind 
plasmas resulting in no dislocation or disruption 
of the streamer. In our study, blobs are thought 
to form along elongating closed arcades below the 
streamer cusp. The closed streamer flux is not in- 
fluenced by the blob formation in any significant 
way. Furthermore, according to our calculations, 
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reconnections occur in the region where the field 
lines are almost parallel to each other, rather than 
in the region where flows/field lines converge most, 
therefore, the idea of converging-flow driven recon- 
nection to form blobs, as proposed by Lapenta & 
Knoll, seems to be inappropriate at least in our 
more realistic streamer configurations. 

Our calculations started from the quadrupo- 
lar magnetic field with a neutral point high in 
the corona, exactly in the form provided by An- 
tiochos et al. (1999) in their break out model 
for CMEs. They investigated the dynamical re- 
sponse of the closed arcades to foot point shearing 
motions. One key point of the breakout model 
lies in the presence of a multipolar pre-eruptive 
field with a magnetic null (or neutral point) high 
in the corona. With the shearing motions, the 
high magnetic null gradually extends into a cur- 
rent sheet structure that manages to survive for a 
certain time in the corona plasmas. Through mag- 
netic forces the current sheet structure serves as 
a major agent confining the increasing free mag- 
netic energy of the sheared central arcades. Erup- 
tions take place when the thickness of the current 
sheet structure eventually decreases to trigger the 
fast break-out reconnections through either phys- 
ical or numerical resistivity. The original break- 
out model did not take into account the effect of 
expanding plasmas on the global magnetic field 
morphology. Such effect is non-negligible espe- 
cially at the heights where the original breakout 
model is concerned and at the null neighborhood 
with very weak magnetic field strength. The con- 
sequence of this effect has been clearly illustrated 
by the calculations presented in this article, espe- 
cially, by comparisons of the upper panels of our 
first figure. We found that the original morphol- 
ogy of the quadrupolar magnetic field is re-shaped 
into a triple-streamer configuration with its neu- 
tral point blown away by the solar wind, indicating 
the critical impact of the solar wind plasmas on 
the coronal field topology that causes the failure 
of the breakout mechanism within the framework 
of our modelling. In future studies, we are going 
to construct a CME model based on the triple- 
streamer solar wind state obtained in this paper 
using the ideal catastrophe mechanism of mag- 
netic flux rope system as the triggering and ma- 
jor energizing mechanism for eruptions (see, e.g., 
Chen et al., 2006, 2007). 
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Fig. 1. — Magnetic field lines superposed by the 
contour map of the two-dimensional distribution 
of the radial velocity field from 1 to 10 R®. The 
field lines are represented by the contours of the 
magnetic flux function with the contour level sep- 
arated by ip = 0.05^0- Panel a) is for the ini- 
tial quadrupolar field and solar wind velocity, b) - 
d) depict the calculated conditions at various in- 
stants. The separatrix associated with the orig- 
inal quadrupolar field (-0 = O.36970o) is shown 
in red curves. The three lines in green, red, 
and yellow colors in the blob region are given by 
= 0.655, 0.66, and 0.665^0, referred to as line A, 
B, and C in the text. To show features in the con- 
tour map of the radial velocity, the magnetic field 
lines are omitted in one of the hemispheres in b) - 
d). [The long-term dynamical evolution (in 2 days 
from t=0 to 4-8 hours) of the triple- streamer solar 
wind blob system is shown in an mpeg animation 
in the electronic edition of the Journal] 
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Fig. 2. — Ten image frames of the side streamer 
part cropped from the corresponding global ver- 
sions of images, spanning from t=19 to t= 23.5 
hours. The frames are separated by 30 minutes 
in time interval with time increasing from left to 
right and top to bottom. The three black arcs 
are centered at the sun with radius being 2.35, 
2.6, and 3.6 Rq. [This figure is available as part 
of the mpeg animation in the electronic edition of 
the Journal] 
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Fig. 3. — The radial components of velocity and 
magnetic field vectors (v r and B r ) together with 
the Alfven speed (va) for the solutions at t = 21.5 
hours left) and 23.5 hours (right). The parameters 
are shown as a function of co-latitude (A = 90° — 
0) for 3 heliocentric distances of 2.35 (solid), 2.6 
(dotted), and 3.6 (dashed) Rq. 
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